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Abstract 

We present an algorithm for constructing the Wilson operator product expan- 
sion (OPE) for perturbative interacting quantum field theory in general Lorentzian 
curved spacetimes, to arbitrary orders in perturbation theory. The remainder in 
this expansion is shown to go to zero at short distances in the sense of expecta- 
tion values in arbitrary Hadamard states. We also establish a number of general 
properties of the OPE coefficients: (a) they only depend (locally and covariantly) 
upon the spacetime metric and coupling constants, (b) they satsify an associativ- 
ity property, (c) they satisfy a renormalization group equation, (d) they satisfy a 
certain microlocal wave front set condition, (e) they possess a "scaling expansion". 
The latter means that each OPE coefficient can be written as a sum of terms, each 
of which is the product of a curvature polynomial at a spacetime point, times a 
Lorentz invariant Minkowski distribution in the tangent space of that point. The 
algorithm is illustrated in an example. 

1 Introduction 

The operator product expansion [HB] (OPE, for short) states that a product of n local 
quantum fields can be expanded at short distances as an asymptotic series, each term of 
which is given by a model dependent coefficient function of the n spacetime arguments, 
times a local field at a nearby reference point 1 y: 

O il {x 1 )O l2 {x 2 )---O ln {xn) ~ ^2c ili2 ...i n k (x 1 ,x 2 , ■ ■ ■ ,x n ) O k {y) . (1) 

k 
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this paper, we will take y = x n but other more symmetric choices are also possible. 
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This expansion has been established in perturbative quantum field theory on Minkowski 
spacetime |40| . and is by now a standard tool, for example in the analysis of quantum 
gauge theories such as QCD. It has also been proven for conformally invariant quantum 
field theories |HHJ E2J EH], and has in fact played a major role in the developement and 
analysis of such theories [HIEE]- Formal mathematical proofs of the OPE have also been 
given within various axiomatic settings (TSJ |2] for quantum field theory on Minkowski 
spacetime. Given the importance of the OPE in flat spacetime, it is of great interest to 
construct a corresponding version of the expansion in curved spacetime. 

In this paper, we present such a construction, within the framework of perturbation 
theory. It is based on the perturbative construction of quantum field theories on globally 
hyperbolic Lorentzian spacetimes which was recently achieved in a series of papers [2H 
122 I2H1 121]; which in turn were based upon key results of Brunetti, Fredenhagen, and 
Kohler EJ. In these papers, the interacting (Heisenberg) quantum fields in curved 
spacetime are constructed as formal power series in the coupling constant (s), that are 
valued in a certain *-algebra of quantum observables. The basic idea how to construct 
the OPE for these interacting fields is as follows: Suppose we have linear functionals 
from the algebra into the complex numbers, that are labelled by an index k enumerating 
the different composite fields, and a reference spacetime point y, and which form a "dual 
basis" to the composite quantum fields in the sense that tyy(Oj(y)) = 5 k j. Now apply the 
functional with label k to the OPE. Then we immediately see that the operator product 
coefficient in front of the k-th term in the sum on the right side of the OPE ought to 
be given precisely by the c-number distribution obtained by applying that functional to 
the product of fields on the left side of the OPE. Below, we will give a perturbative 
construction of such a dual basis of functionals in the context of a scalar, renormalizable 
field theory model in any 4 dimensional Lorentzian spacetime. In this way, we will obtain 
the desired perturbative formula for the OPE coefficients, C. 

While this construction gives a conceptually clean derivation of a perturbative expres- 
sion for the coefficients C, it does not yet show that the remainder of the OPE expansion 
defined in this way actually goes to zero (and in what sense) when the points are scaled 
together. To analyze this question, we apply any Hadamard state to the remainder in 
the OPE. We then show that the resulting c-number distribution in n points goes to zero 
in the sense of distributions when the points are scaled together in a arbitrary fashion. 
Thus, the OPE holds in the sense of an asymptotic expansion of expectation values, to 
arbitrary order in perturbation theory, and for any Hadamard state. The proof of this 
statement mainly relies on the known scaling properties of the various terms in the pertur- 
bative series for the interacting fields [2H 1221 E] • However, these properties by themselves 
are not sufficient, for the following reason. The perturbative formulae for the interacting 
fields at k-th. order involve an integration over k "interaction points" in the domain of the 
interaction (which might be the entire spacetime). It turns out that we can control the 
contributions from these integrations in the OPE if we split the interaction domain into 
a region "close" to the points xi, . . . ,x n , and a region "far away." However, the points 
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xx,...,x n themselves are supposed to be scaled, i.e., they move, so the split of the inte- 
gration domain has to be constantly adapted. We achieve this by dividing the interaction 
domain into slices of thickness 2~- ? centered about y, where j = 1, 2, 3, ... . We find that 
the contributions from these slices can be controlled individually, and then be summed, 
if the interaction is renormalizable. 

Finally, we will derive the following important general properties about the OPE 
coefficients C. (a) They have a local and covariant dependence upon the spacetime 
metric, (b) they satisfy an associativity property, (c) they satsify a renormalization group 
equation, (d) they satisfy a certain microlocal condition on their so-called wave front set, 
and (e) they can be expanded in a "scaling expansion" . Let us explain these properties. 

The local covariance property (a) of the OPE coefficients states that if (M, g) and 
(M', g') are globally hyperbolic spacetimes with corresponding OPE coefficients C resp. 
C, and if / : M — > M' is a causality and orientation preserving isometric embedding, 
then f*C is equivalent to C at short distances. Thus, in this sense, the OPE coefficients 
are local functionals of the metric (and the coupling constants), and in particular do not 
depend upon the large scale structure of spacetime, such as the topology of M. This 
property can be understood from the fact that, as shown in j23El, the interacting fields 
may be constructed in a local and covariant fashion. More precisely, whenever we have a 
causality preserving isometric embedding /, there exists a linear map aj from the quantum 
field algebra associated with (M, g) into the algebra associated with (M', g') preserving the 
algebraic relations. Furthermore, the fields are local and covariant in the sense |2UE] that 
the image of an interacting field on (M,g) via a/ corresponds precisely to the definition 
of that field on (M',g'). Since the OPE may be interpreted as an asymptotic algebraic 
relation, and since the local and covariance property as just stated means that algebraic 
relations only depend on the metric locally and transform covariantly under spacetime 
embedding, it is natural to expect that also the OPE coefficients depend locally and 
covariantly upon the metric (and the couplings). This is indeed what we shall prove. 

The associativity property (b) arises when one studies the different ways in which a 
configuration of n points can approach the diagonal in M n . For example, in the case n — 3, 
we may consider a situation in which all three points approach each other at the same 
rate, or we may alternatively consider a situation in which two points approach each other 
faster than the third one. The possible ways in which n points can approach each other 
may be described by corresponding merger trees T which characterize the subsequent 
mergers [T51 IT]. The associativity property states that scaling together the points in an 
operator product according to a given tree is equivalent to performing subsequently the 
OPE in the hierarchical order represented by the tree T. We will argue — relying mainly on 
a general theorem of j2Hj — that this type of "short distance factorization property" indeed 
holds in perturbation theory. The associativity may again be understood intuitively from 
the fact that the OPE coefficients are in some sense the structure "constants" of the 
abstract associative *-algebra of which the interacting fields are elements. 

The renormalization property (c) of the OPE coefficients arises from the fact that the 
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algebras of quantum fields satisfy a similar property [23] ■ Namely, if we rescale the metric 
by a constant conformal factor A 2 , then this is equivalent (in the sense of giving rise to 
isomorphic algebras) to redefining the field generators in a particular way, and at the 
same time letting the coupling constants of the interaction (which enter the structure of 
the algebra) flow in a particular way dictated by the "renormalization group flow" of the 
theory [2Sj- Again, since the OPE coefficients are in a sense the structure constants of 
the algebra, they can be expected to have a corresponding property. 

The microlocal property, (d), is a property describing the nature of the singularities 
in an OPE coefficient C. It states that the wave front set [2Z| WF(C) has a characteristic 
form that encodes the positivity of energy momentum in the tangent space. Our condition 
found for WF(C) is similar in nature to the so-called "microlocal spectrum condition" 
ES] for the wave front set of correlation functions of linear field theories in curved space. 
However, our condition differs qualitatively from that proposed in |33] in that the 
interactions may affect the form of the wave front set in our case. 

The scaling expansion (e) states that, if we scale n points xi, X2, ■ ■ ■ , x n in M together 
according to a merger tree, then a given OPE coefficient C can be approximated to any 
desired precision by a finite sum of terms each of which has the form of a polynomial in 
the mass and curvature tensors at point y = x n , times a Lorentz invariant Minkowski dis- 
tribution in the Riemannian normal coordinates of X\, . . . , x„_i relative to y = x n . These 
Minkowski distributions have, to each order in perturbation theory, a simple homogeneous 
scaling behavior modified by polynomials in the logarithm. They may be extracted from 
the given OPE coefficients by taking a certain "Mellin moment", which is an operation 
defined by first taking the Mellin-transform of a function and then extracting certain 
residue. 

The general properties just described (except (e)) are postulated axiomatically in the 
fourthcoming paper [22], arid so our present analysis may be viewed as a confirmation of 
these axioms in perturbation theory. 

Our plan for this paper is as follows. In section 2, we first recall the general strategy for 
obtaining the perturbation series for the interacting fields [2U 122111]. The OPE is derived 
in section 3, and its general properties are derived in section 4. An example illustrating 
our algorithm for computing the OPE coefficients is presented in section 5. 

As indicated, for simplicity and concreteness, we only consider a single hermitian, 
scalar field with renormalizable interaction (in 4 spacetime dimensions). While the re- 
striction to a renormalizable interaction seems to be essential, we expect that our algo- 
rithm will work for other types of fields with higher spin and renormalizable interaction, 
in other dimensions. However, the analysis of the OPE in the physically interesting case 
of Yang-Mills theories would first require an understanding of the renormalization of such 
theories in curved spacetime, which is considerably complicated by the issue of gauge 
invariance. This is at present an open problem. 
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2 Perturbation theory 



A single hermitian scalar field in 4 dimensions is described classically by the action 



S 



dfi . (2) 



+ (m 2 + aR)<p 2 + 2 V 

' M I 

Here dfi = ^Jg dx° A ■ ■ • A dx 3 , m 2 , a G R, the quantities «j G R are coupling parameters 
parametrizing the strength of the self-interaction of the field, and throughout this paper 
Oi are polynomials in the field (f) and its covariant derivatives, as well as possibly the 
Riemann tensor and its derivatives. In the above action, they encode the nature of the 
self-interaction. Later, we will assume that the interaction is renormalizable, but for the 
moment no such assumption need to be made. 

The perturabative construction of the quantum field theory associated with this action 
has been performed in a series of papers |2"H l2"2"| E3J Ej • These constructions consist of the 
following steps. First, one defines an abstract *-algebra f2J HI] ^{M,g) containing the 
quantized field 0, together with its Wick powers Oi, for the corresponding linear theory, 
which classically corresponds to dropping the self-interaction term 



I — 2_] KiOi dfx = Cd\± 

J M „• J M 



(3) 



from the above action S. From these quantities, one then constructs the corresponding 
interacting quantum fields (smeared with a compactly supported testfunction) as formal 
power series in free field quantities via the Bogoliubov formula |T!H I29j. 

i n / \ 

Oi(h)j = J2^ R n[LA h d ^ J0 "j » h e C ™( M ) » ( 4 ) 

n>0 

or more formally without smearing, 

•n 

O i (x) I = y2-R n (O l (x),I® n ), (5) 

n>0 

where the quantities R n are the so-called retarded products, which are multi-linear maps 
on the space of local classical action functionals of the form Q, taking values in the 
underlying free field algebra JF(M, g) . In order to describe these constructions in more 
detail, it is first necessary to recall at some length the definition and key features of the 
linear field algebra JF(M, g) and its quantum states, as well as nature of the retarded 
products R n . 

The definition of T{M,g) can be stated in different equivalent ways [2~3 | WI \ HI], 
and we now give a definition that is most suited to our purpose. Let u)2 G T>'(M x M) be a 



5 



bidistribution of Hadamard type, meaning (a) that uj 2 is a bisolution to the Klein-Gordon 
equation □ — m 2 — aR in each entry, (b) that the anti-symmetric part of u 2 is given by 

lu 2 {xi, x 2 ) - u 2 (x 2 , si) = i(A A (x 1 , x 2 ) - A R (xx,x 2 )) = iA(xx,x 2 ) , (6) 

where Aa/r are the unique advanced and retarded propagators of the Klein-Gordon equa- 
tion ^7j, and (c) that it has a wave front set of Hadamard type |35j : 

WF(u 2 ) = {(>!, h; x 2 , k 2 ) eT*(MxM)\ 0; k x = p, k 2 = -p,p G V;} . (7) 

Here, it is understood that the set can only contain those x\ and x 2 that can be joined 
by a null geodesic, 7, and that p = p^dx^ is a parallel co-vector field tangent to that 
null geodesic, meaning that V^p = 0, and V± is the dual of the future/past lightcone. 
The desired Wick-polynomial algebra T{M, g) is now generated by an identity 1, and the 
following symbols F: 

P n 

F= u n )(x u ...,x n ) A d/ifo) , (8) 

where symmetric, compactly supported distribution on M n subject to the wave 

front condition ^2] 

WFK)n[(y;) n u(K!) n ]=0. (9) 

The relations in the *-algebra T{M, g) are as follows: The *-operation is defined by letting 
F* be given by the same expression as F, but with u n replaced by its complex conjugate, 
and the product is defined by 

• -ui ■ ■ ■ ; Xyi) . (j) .(j {x n -\-\ 1 • • • , ^ri+rn) 

k K J ' K >' ' Pl ,..., Pk eP i 

where P is the set of all pairs pi G {1, . . . , n} x {n + 1, . . . , n + m}. This formula is identical 
in nature to the standard Wick theorem for normal ordered quantities (relative to a Gaus- 
sian state with 2-point function uj 2 ). The wave front conditions on u n and uj 2 are needed 
in order guarantee that the product between the corresponding integrated quantities as 
in (jHJ) exists, because the latter involves the pointwise products of distributions |2Z] • The 
definition of the algebra T{M, g) superficially seems to depend on the particular choice 
of o>2, but this is in fact not so: A change of u 2 merely corresponds to a relabeling of the 
generators, and does not change the definition of T{M,g) as an abstract algebra fJT]. 

The relation between the abstract quantization of the linear Klein-Gordon field just 
described and more familiar ones is as follows. If / is a smooth test function, then the 
generator (jHJ) with n = 1 and u\ = f can be identified with the smeared field </>(/). Indeed, 
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using the Wick formula (JlOjh and the antisymmetric part of u>2, eq. (JHJ), one easily derives 
the relation [<fi(fi), 4>(f2)] = /z)!, which is the standard commutation relation for 

a linear scalar field in curved spacetime. The higher order generators (JHJ) with u n given 
by the n-fold tensor product /® n = / <g> • • • <g> / of a smooth testfunction correspond to a 
smeared normal ordered product : cj)® n : u {f® n ), formally related to the field itself by 



5 n f 1 



■ (11) 



As it stands, the smeared version of the Klein-Gordon equation </>((□ — m 2 — aR)f) = is 
not an algebraic relation. However, this relation could easily be incorporated by factoring 
J-{M, g) by the 2-sided *-ideal J{M, g) consisting of all elements F of the form (JHJ) with 
u n a distribution in the class (JHJ) which is in the image of this distribution class under the 
Klein-Gordon operator, such as (□ — m 2 — aR)f in the simplest case. The purpose of this 
paper will be to establish an OPE, and it is technically convenient for this purpose not 
to factor by the ideal. However, after the OPE has been constructed there is absolutely 
no problem to factor by this ideal, because it is clear that the OPE will continue to hold 
on the factor algebra. 

Quantum states in the algebraic framework are linear expectation functionals $ : 
J-{M,g) — > C that are normalized, meaning $(1) = 1, and of positive type, meaning 
$>(F*F) > 0. Of particular importance are the so-called Hadamard states on J-"(M,g). 
Those states are defined by the fact that their 2-point function $2(^1, x 2 ) = §(4>(. x i) ( f ) { x 2)) 
satisfies properties (a), (b), and (c) listed above, and that their truncated n-point functions 
of the field </> for n 7^ 2 are smooth solutions to the Klein-Gordon equation. The key 
consequence of the Hadamard requirement which we shall need later is J20j that $(: (f)® n : w 
(x\, . . . ,x n )) is smooth. Note that by definition, the n-point functions of a Hadamard 
state satisfy the Klein-Gordon equation. Consequently, they vanish on the ideal J7"(M, g) 
generated by the Klein-Gordon equation and hence induce states on the factor algebra. 

Later, we want to define an operator product expansion, and for this we will need a 
notion of what it means for algebra elements to be "close" to each other. For this we 
now introduce a topology on T{M,g). There are various ways to do that. A particular 
topology was introduced in We prefer here to work with a different (weaker) topology, 
defined by the collection of all linear functionals on JF(M, g) with the property that 
$(: 0® n : w (xi, . . . ,x n )) is smooth. This set includes the Hadamard states as defined 
above, and we shall, by abuse of notation, sometimes refer to such $ as "Hadamard" as 
well. We then introduce a set of seminorms f\f$(F) = labelled by these functionals 

$. We say that a sequence {-F/vj-jvgN of algebra elements tends to zero if for each $, and 
each e > 0, there is an N such that J\f®(F N ) < £ for all iV > N . 

An important feature of the algebra T{M, g) is that it has a local and covariant 
dependence upon the spacetime. More precisely, if / : M — > M' is a causality and 
orientation preserving isometric embedding of a globally hyperbolic spacetime (M, g) into 
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another such spacetime (M', g'), then there exists a continuous, injective *-homomorphism 



a f '■ F(M, g) — ► T{M' , g') . (12) 

This embedding is most simply described in terms of its action on a smeared field <p(h), 
where h is a test function on M. If hi = f*h is the corresponding pushed forward test 
function on M', then we define af[<fi(h)] = <f)(h'). Furthermore, the action of aj on an 
arbitrary element in T{M, g) may then be defined by continuity, because the subalgebra 
generated by expressions of the form <p(h) is dense in J^(M, g). The action of o.f on 
the smeared field (f>{h) is characteristic for so-called "local covariant fields". Namely, an 
algebra valued distribution Oi : Cg°(M) — > J-(M, g), h i— > Oi(h) that is defined for all 
spacetimes (M, g) is called a (scalar) local and covariant field if 

a f [O i (h)\ = O i {ti) h' = f*h, (13) 

whenever / is an orientation and causality preserving isometric embedding. Local covari- 
ant fields of tensor type are defined in the same way, except that the testfunction h is now 
a section in the (dual of the) vector bundle Ej corresponding to the tensor type. Thus, 
the field is (by definition) a local and covariant field. On the other hand, the normal 
ordered n-th Wick power of a field defined by putting u n (xi, . . . , x n ) = f(x n )S(xx, . . . , x n ) 
in eq. (jHJ) is not a local and covariant field, because it implicitly depends on the choice 
of the 2-point function co> 2 , which is not a local and covariant quantity [21]. The possible 
definitions of Wick powers giving rise to local and covariant fields (satisfying also various 
other natural conditions) were classified in [21] . It turns out that the definition of a given 
classical expression 

d = V ai 0V a2 • • • V a "0, % = {a a , a 2 , . . . , a n } e 1 = Z| (14) 

n 

as a local, covariant field in F{M, g) is not unique, but contains certain ambiguities. As 
proven in [21], these ambiguities correspond to the possibility of adding to a given field 
lower order Wick power powers times certain polynomials of the Riemann tensor and its 
derivatives V( Ml • • • V^)^?^^^^ of the same dimension as Oi. Here, the dimension of a 
Wick power is the map [•] : X — > N from the index set labelling the various fields, into the 
natural numbers defined by 

n 

[i]=n + J2 a i- ( 15 ) 
i=i 

One definition which is local and covariant (and satisfies also the other natural conditions 
given in |21j ) is the following "local normal ordering prescription". It is based upon the 
use of the local Hadamard parametrix H, which is the bidistribution defined on a convex 
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normal neighborhood of the diagonal {(x,x); 16 M} of M x M by 2 



\n>0 

In this expression, a(x,y) is the signed squared geodesic distance between two points, 
we have defined t(x,y) = r(x) — r(y), where r is a time-function, and the v n (x,y) are 
smooth symmetric functions that are determined by requiring that if be a parametrix, 
i.e., a solution to the Klein-Gordon equation □ — m 2 — aR in each entry modulo a smooth 
remainder. Explicitly v = D 1 ^ 2 /2tt 2 is given in terms of the VanVleck determinant D, 
defined by 

4 \J\ 1 ; 

Here, J(x,y) G T*M <g> T*M is the bitensor of parallel transport, (V <8> V)a(x, y) G 
T*M ® T*M is the is the bitensor obtained by taking the gradient of a(x, y) in both x 
and y, and we are defining a biscalar \B(x,y)\ by 

4 

\B(x, y) | dnix) ® dn(y) = f\ B(x, y) (18) 

from any bitensor B(x,y). The smooth functions v n (x,y) are iteratively defined by the 
transport equations [TO] 

2(W) V^ n + [( W) V„ In D + An}v n = -2(D - m 2 - aR)v n ^ , (19) 

where the derivatives act on x. These functions are symmetric in x and y [SUEI, and 
their germs at the diagonal are locally and covariantly defined in terms of the metric. 
Where it is well-defined, H has a wave front set of Hadamard type (|7J). Next, fix a convex 
normal neighborhood of the diagonal in M n , and in that neighborhood define locally 
normal ordered products : 4>® n : H (xi, . . . ,x n ) by the same formula as (JUJ), but with u; 2 
in that formula replaced by H, 



■ (20) 



As the expressions pijl . their expectation value in any Hadamard state is smooth. Fol- 
lowing [21] we define the local covariant n-th Wick power of the field as the distribution 
valued in T{M,g) given by 

<f> n (x) = lim : <P® n : H (x(e)) x(e) = (exp^i), . . . , exp^e^n)) , (21) 

£^0 



2 The infinite sum is to be understood in the sense of an asymptotic expansion. 
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where £j denote the Riemannian normal coordinates (identified with a vector in R 4 ) of 
the point Xi relative to x. More generally, fields containing derivatives are defined by 
first acting with the derivatives on the appropriate tensor factor in : (p® n : H before taking 
the above "coincidence limit". That is, if i ) G X denotes a collection of 

natural numbers, then the corresponding local covariant field Oi(x) is defined by applying 
the partial derivative operator ■ ■ -<9g™ prior to the coincidence limit. Note that this 
definition is covariant, because partial derivatives with respect to Riemannian normal 
coordinates at x may be expressed in terms of curvature tensors at x and covariant 
derivatives V. 

Having described the algebra T[M } g) of field observables in the linear field theory 
associated with the action 5* without the interaction terms, we now turn to the interacting 
theory. For this it is technically convenient at an intermediate step to assume that the 
couplings Ki in the action S are not constants, but actually smooth functions of compact 
support in M, which we assume are locally constant, 

Ki(x) = Kix(x) , (22) 

where \ G C£°(M), and x{ x ) = 1 in an open set in M with compact closure. The 
cutoff functions x serve as an infra red cutoff and are removed at a later stage. With 
the introduction of cutoff functions understood, the interacting fields are defined by the 
Bogoliubov formula © in terms of the retarded products R n . Each retarded product is 
a continuous, bilinear map 



Rn-Sx (9)5 \^F{M,g) (J ,/i® ■■■®I n ) »Rn(I ,h® ■■■®In), (23) 




from the tensor powers of the space 5 of all classical action functionals Ij = j Cj dfj, that 
are local and polynomial in the field <fi and whose couplings are compactly supported 
functions on M. The map is taking values in the algebra JF(M, g) associated with the 
linear field theory, and is symmetric in I\, . . . , I n . Note that the power series expression (jSJ) 
for the interacting fields is only a formal series, and no statement is made about its 
convergence. Since each term in these series is an element in !F(M,g), the interacting 
fields are elements of the algebra V £g> T{M, g), where 

V = C[K 1 ,K 2 ,...] = \^2a ai ... ah K^---K^; a ai ... ak ec[ (24) 

U t >o J 

is the corresponding ring of formal series. All operations, such as multiplication in this 
ring and in the algebra V <g> T{M, g), are defined by simply formally multiplying out the 
corresponding formal series term by term. Furthermore, this algebra inherits a natural 
topology from !F(M,g): A formal power series converges to another formal power series 
in the algebra if each coefficient does. 
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We want the zeroth order contribution R (Oj(x)) in the Bogoluibov formula (0) to be 
given by our definition of Wick powers Oi(x) in the linear field theory 3 . Thus, we define 
Ro(Oi(x)) to be equal to the locally normal ordered field Oi(x) given in eq. (J2~T| . The 
terms with n > 1 in the formal series (JHJ) represent the perturbative corrections coming 
from the interaction, I. They involve the higher, non-trivial, retarded products. The con- 
struction of these retarded products can be reduced 4 to the construction of the so-called 
"time-ordered products," because there exists a well-known formula of combinatorial na- 
ture relating these two quantities, see e.g. the appendix of [Hj. The construction of the 
time-ordered products in turn has been given in [2D E2J Elj , which is based on work of j3j. 
The strategy in these papers is to first write down a number of functional relations for the 
time ordered products that are motivated by corresponding properties of the interacting 
fields defined by the Bogoliubov formula. These properties then dictate to a large extent 
the construction of the time-ordered (and hence the retarded) products. Since there is a 
combinatorial formula relating the time ordered products to the retarded products, these 
relations can be equivalently be stated in terms of the retarded products. The relevant 
relations 5 for this paper are as follows: 

(rl) Causality: Let F,G,Si G S. Suppose that there is a Cauchy surface such that 
suppG is in its future and suppF in its past. Then 



(r3) Expansion: There exist local covariant c-number distributions r near the total 

3 Note that the argument of the retarded product Ro(0i(a;)) is a classical action (or density), while 
the corresponding Wick power Oi(x) is a distribution valued in the quantum algebra !F(M, g). We should 
strictly speaking distinguish these quantities by introducing a new notation for the Wick power, but we 
shall not do this for simplicity. 

4 Alternatively, it should also be possible to construct the retarded products directly along the lines 
of [1211131 HI) , by suitably generalizing the arguments of that paper from Minkowski spacetime to curved 
spacetime. 

5 A complete list may be found in |24j . 




(25) 



(r2) GLZ factorization formula (Hj: 




(26) 
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diagonal {(x, x, . . . , 



x)] x G M} such that (with lj < ij) 



R n (Oi (x), O ix {y\) ■ ■ ■ O in (y n )) 



ioh— in 

!o,il,...,Zn 



(x,yi,...,y n ) : Oi (x)O h {yi) . . .Oi n {y n ) : 



ii 



(27) 



(r4) Scaling degree: The distributions r have the scaling degree 



sd(rj°t.t) = £[*] " M 



(28) 



at the total diagonal {(x, x, . . . , x); x G M} C M n+1 . Here, the scaling degree of a 
distribution u G T>'(X) at a submanifold 7 C I is defined as follows [HUE]- Let 
S £ : X — ► X be an injective, smooth map defined on an open neighborhood X 
of Y with the properties (a) that S e \y= idy, and (b) that for all y G Y, the map 
(DS s )(y) : T y X T y X is the identity on T y Y and scales vectors by e > on a 
complementary subspace C y C T y X of T^y. Then u has scaling degree sd(tt) at Y 
if lim £ _ >0 ^ D ' u ° Se = f° r ai l D > sd(tt) in the sense of X>'(Xq). The definition is 
independent of the precise choice of S E . 

(r5) Locality and covariance: Let / : (M,g) — > (M',g') be a causality and orientation 
preserving isometric embedding. Then the retarded products satisfy 



where /* denotes the natural push-forward of a local action functional on M to the 
corresponding action functional on M'. 

(r6) Microlocal condition: The distributions r in the expansion ()27|) have the following 

wave front set: 



WF(r) C < (x, k; y%, 1%; . . . ; y m , l m ) G T*M m+ ; there is graph in C/i )n such that 



The valence of the vertices yt in the graph is restricted to be less or equal than the 
maximum power of occurring an the operators Oi l , . . . , Oi n in eq. ()27j1 . 




(29) 




e:t{e)=yi 



Pe 



(30) 
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In the formulation of the last condition, we are using a graph theoretical notation [3], 
which will be useful later as well. Most generally, we consider the set G n ,m of embedded, 
oriented graphs in the spacetime M with n + m vertices. Each such graph has n so-called 
"external vertices", G M, and m so-called "internal" or "interaction vertices" 

Hi, ■ ■ ■ , Vm £ M. These vertices are of arbitrary valence, and are joined by edges, e, which 
are null-geodesic curves 7 e : (0, 1) — > M. It is assumed that an abstract ordering 6 < of 
the vertices is defined, and that the ordering among the external vertices is x\ < ■ ■ ■ < x n , 
while the ordering of the remaining interaction vertices is unconstrained. If e is an edge 
joining two vertices, then s(e) (the source) and t(e) (the target) are the two vertices 7 e (0) 
and 7e(l), where the curve is oriented in such a way that it starts at the smaller vertex 
relative to the fixed vertex ordering. Each edge carries a future directed, tangent parallel 
covector field, p e , meaning that V^ e p e = 0, and p e G dVT. 

Similar to the case of local covariant Wick products (the case n — 0), the above func- 
tional relations (together with other functional relations described in detail in [2D 122]) do 
not uniquely fix the retarded products: There remains a number of real constants at each 
order n which parametrize the set of possible definitions of R„ that are compatible with 
(rl)-(r6). These correspond to the usual "renormalization ambiguities" in perturbative 
quantum field theory, see [221 123 121] for details. 

One finally needs to remove the dependence of the interacting fields on the arbitrary 
cutoff function \. For this, one investigates how the interacting field changes when the 
cutoff function is varied. Assume that Xi an d X2 are two different cutoff functions, both 
of which are equal to 1 in an open globally hyperbolic neighborhood U M . Let I\ and 
I2 be the corresponding interactions. Note that, as a classical functional, the difference 
I\ — I2 is supported in a compact region, and vanishes in the neighborhood U where 
the cutoff functions coincide. The key fact |3], which follows from the above functional 
relations (rl) and (r2), is now that there exists a unitary operator V G V <8> J-(M,g), 
depending upon J 1( J 2 , with the property that 

Oi(x) h = V Oi(x) l2 V*, for all x G U, i G J. (31) 

This relation may be interpreted as saying that the algebraic relations between the inter- 
acting fields within the region U where the cutoff function is constant do not depend on 
how the cutoff function is chosen outside this region, and this observation may be used to 
construct an abstract interacting field algebra associated with the entire spacetime (M, g) 
that is independent of the choice of cutoff function [3]. However, in the present context, 
we are actually only interested in a small patch U of spacetime where we want to consider 
the OPE. Therefore, it will be more convenient for us to simply fix an arbitrary cutoff 
function that is equal to 1 in the patch U of interest. Since the OPE concerns only local 
algebraic relations, it is clear from (|31j) that it should not matter what cutoff function we 
choose, and this will formally be shown below in item 2) of section 4. 

6 The ordering is not assumed to be related to the causal structure of the manifold at this stage. 
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3 Operator product expansion 



We will now show that the interacting fields Oi(x)i described in the previous section obey 
an operator product expansion, 

Oi l (x 1 )iO ia (x i )i---Oi n (x n ) I ~ ^2 C ili2 ... in k (x 1 ,x 2 ,...,x n ) I O k (x n ) I , (32) 

[Jfe]<A 

where Cj are certain distributinal coefficients depending upon the interaction, J, which 
are to be determined, and where [k] is the standard dimension function defined above. 
We mean by the above expression that, as the points x\, x 2 , ■ ■ ■ , x n approach each other, 
the algebra product of the interacting fields on the left side can be approximated, to 
the desired precision determined by A, by the right side, in the topology on the algebra 
V®F{M,g). 

To make this statement precise, we must, however, take into account that both sides 
of the OPE are actually distributional, and that a configuration ) e M n 

of n mutually distinct points on a manifold may "merge" in qualitatively different man- 
ners when n > 2, because the points may approach each other at "different rates". The 
appropriate mathematical framework to formalize in a precise manner the possibility of 
configuration of points to approach each other at different rates is provided by a con- 
struction referred to as the "compactification of configuration space," due to Fulton and 
MacPherson |18j . and Axelrod and Singer PQ. Let 

M n = {x=(xi,X2,...,x n )eM n ;x i ^x j } (33) 
= {x G Map({l, . . . , n}, M), x(i) — X{] x inject ive} (34) 

be the configuration space, i.e., the space of all configurations of n mutually distinct points 
in M. The union of partial diagonals 

dM n = |J A s C M n (35) 

Sc{l,...,n} 

where a partial diagonal is defined by 

A s = {x G Map({l, . . . , n}, M); x \ S = constant} (36) 

is the boundary of the configuration space Mq. Configurations of points where some points 
come close to each other are in some sense close to this boundary. The Fulton-MacPherson 
compactification M[n] is obtained by attaching a different boundary, dM[n], to Mq, which 
in addition incorporates the various directions in which dM£ can be approached. This 
boundary may be characterized as the collection of endpoints of certain curves x(e), in 
M[n], which are in Mq for e > 0, and which end on dM[n] at e — 0. These curves are 
labeled by trees T that characterize subsequent mergings of the points in the configuration 
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as £ — > 0. A convenient way to describe a tree T (or more generally, the disjoint union of 
trees, a "forest") is by a nested set T = {Si, . . . , S^} of subsets Si C {1, . . . , n}. "Nested" 
means that two sets are either disjoint, or one is a proper subset of the other. We agree 
that the sets {1}, . . . , {n} are always contained in the tree (or forest). Each set Si in T 
represents a node of a tree, i.e., the set of vertices Vert(T) is given by the sets in T, and 
Si C Sj means that the node corresponding to Si can be reached by moving downward 
from the node represented by Sj. The root(s) of the tree(s) correspond to the maximal 
elements, i.e., the sets that are not subsets of any other set. If the set {1, . . . ,n} G T, 
then there is in fact only one tree, while if there are several maximal elements, then 
there are several trees in the forest, each maximal element corresponding to the root 
of the respective tree. The leaves correspond to the sets {1}, . . . , {n}, i.e., the minimal 
elements. 

The desired curves x(e) tending to the boundary of M[n] are associated with trees 
and are constructed as follows. With the root(s) of the tree(s), we associate a point 
Xi G M, where i is a label that runs through the different maximal elements, while with 
each edge e G Edge(T) of a given tree (a line joining two nodes), we associate a vector 
v e G T Xt M, where is associated with the root of the tree that e belongs to. To describe 
the definition of this vector it is convenient to identify an edge e G Edge(T) with the pair 
e = (S, S') of nodes that it connects, i.e., an edge defines a relation in T x T. If S C 5", 
then we write S' = t(e) for target, and S = s(e), for the source. We then set 

Ve = £m(t(e)) - fm(s(e)) TTl(S) = maxjij i E S} , (37) 

where £j denotes the Riemannian normal coordinates (identified with a vector in R 4 via 
a choice of orthonormal tetrad at the corresponding root Xj). We define the desired curve 
x(e) by 

Xj (e) = Y, v e£ deptHt{e)) ■ (38) 

where pj is the unique path connecting the leaf j with the corresponding root, where 
v e is given in terms of x by eq. (J37J), and where depth(S) is the number of edges that 
connect the node S G T with the root. The following figure illustrates this definition in 
an example: 
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For each fixed tree T, and each fixed e, the above curve defines a map 

^r(e) : M n -> M n , x i-> £(e) (39) 

flowing the point x = x(l) to the point x(e). For e = 0, the image of this map may be 
viewed as a portion of the boundary 9M[n] corresponding to the tree. The roots of T 
correspond to the particular diagonal; in particular, if there is only one tree in T (as we 
shall assume from now on) then the configuration x(e) converges to the total diagonal in 
M n . It may be checked that the maps ijirie) satisfy the composition law 

M^°M^) = M^')- ( 4 °) 

Using the maps ^>r{s) we can define an asymptotic equivalence relation ~s,t for distribu- 
tions on M n . Consider distributions u±, u 2 defined on M n . For a given tree T and 5 > 0, 
we declare the equivalence relation ~r,<5 by 

u\ ~r,5 u 2 lim e~ s (u t - u 2 ) o V'r(e) = , (41) 

in the sense of distributions on M n , where we view as a map M n — > M n that is 

parametrized by e > 0. 

Having defined the equivalence relation we can now state precisely our notion of 
an OPE. Namely, we require that, for each 5 > 0, each given set of operators, and each 
tree T, there exists a A so that the OPE holds in the sense of ~5,r- The only issue that we 
have not yet been quite precise about is that the OPE is not a relation between c-number 
distributions, but instead distributions valued in the topological algebra P®JF(M, g). This 
difficulty is simply dealt with by requiring convergence in the equivalence relation (|4ip 
(now for algebra valued objects) with respect to the topology in the algebra. Thus, we 
define the precise sense in which the OPE is supposed to hold to be that for each tree T 
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with one root, and each 5, there exists a A G R such that (J52j) holds in the sense of ~5 5 r 
as a relation between the corresponding algebra valued distributions. 

We now come to the actual construction of the operator product coefficients in per- 
turbation theory. As also described in j23- and as originally suggested by Bostelmann j2j 
and Fredenhagen and Hertel ^H] in the context of algebraic quantum field theory on 
Minkowski spacetime, it is convenient to think of the operator product coefficients as 
arising via certain "standard functionals" 

*Af,x( ■)r-V® ?{M, g)^V®E i \ x . (42) 

These functionals depend upon the given spacetime (M, g), the index label i Gl describ- 
ing a composite field, a point x G M, and the interaction, as indicated by the subscript 
"J" . The functionals take values in the fiber over x in the vector bundle Ej (viewed as a 
P-module) associated with the tensor character of the field 0j. In our constructions be- 
low, the functionals are in fact only defined on the subalgebra V ® ^{U, g) corresponding 
to a convex normal neighborhood U C M. However, since all of our considerations are 
entirely local, we may assume without loss of generality and to save writing that U = M. 

The OPE coefficients are supposed to be given in terms of the above standard func- 
tionals by 

C lll2 ...J{x u x 2 , . . . ,z n )j = ¥ Xn ((!\(xi)/e> i2 (x2)j ■ • • CUO/)j • (43) 

We will construct the OPE coefficients in perturbation theory by presenting a suitable 
set of such standard functionals. We are going to choose these standard functionals as a 
"dual basis" to the interacting fields, in the sense that we wish them to satisfy 

^[.{Pj{x)i)i = 6 { j id El for all x G M and [i], [j] < A. (44) 

This ansatz is motivated by the following simple consideration. Let us assume that an 
OPE exists. Let us fix a A > 0, carry the OPE out until [k] < A, and apply the functionals 
\fr 3 x to it, where [j] < A. Using (j44j) . we immediately find that the coefficients in the 
OPE must be given by (|4H|) . up to a remainder term coming from the remainder in the 
OPE. But this remainder is by assumption small for asymptotically short distances, in 
the sense of the above equivalence relation, provided we make A sufficiently large. It can 
therefore be ignored. 

Thus we have argued that if an OPE exists in the sense above, and if standard func- 
tionals satisfying (|4*4"|) have been defined, then the OPE coefficients C ought to be given 
by AH)) . Consequently, our first step will be to define the standard functionals as formal 
power series in the coupling constants so that eq. will be satisfied to arbitrary 
orders in perturbation theory. To zeroth order in perturbation theory, such standard 
functionals are defined as follows (see also [2S1)- Recall that a general algebra element 
F G !F(M,g) can be written as in eq. (JHJ) in terms of normal ordered generators If 
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we are interested only in elements F so that the corresponding u n in (JSJ) are suppored 
sufficiently close to the diagonal in M n (as we will always assume in the following), then 
we may rewrite F in terms of the locally normal ordered generators : <fi® n : H • • • , x n ) 
given in eq. ()2()j) instead of the normal ordered generators : <p® n : w (x\, . . . , x n ). The 
action of the zeroth order standard functionals is then declared by 

% (: ^ m : H (*!, . . . , x m )) = 5m > n ef ai " " " C" i = (ai,...,a n ), (45) 

Oi - ■ ■ ■ o n . 

and extended to all of T{M, g) by linearity. Here, £j are the Riemannian normal coor- 
dinates of Xi relative to x, identified with vectors in T X M. These functionals satisfy the 
analog of eq. (jBj) for the linear fields defined above. Since the interacting fields Oi(x)j 
are given by formal power series whose zeroth order is the linear field expression (see the 
Bogoliubov formula ©), it follows that the action of the linear field functionals on an 
interacting field is of the form 

*t{O j (x) I ) = 5 l j ii 1k +A i j (x), (46) 

where A l j is the endomorphism in End(Ej,Ej) that arises from the higher perturbative 
contributions to the interacting field, see (J5J), and is given by 

•n 

n>l 

Consequently, using the standard geometric series for the inverse of a linear operator of 
the form 1 + L and writing out explicitly the above formula for A l j(x), we find that the 
functional defined by the following series is a solution to the equation (j44j): 



miH \-m k 



rav. - ■ ■ rau\ 
fc=0 mi>l 

V (RmAOn; / 0mi )) (R m2 (O j2 ; /® m2 )) • • • V k (R mk (0 3k+1 ; I® m *)) ^ (F) . (48) 

Here, m = mi is the perturbation order of an individual term, and the sums over j\ 
are carried out to order < A. Thus, for each fixed m, the sum over k has only a finite 
number of terms, and the resulting expression is a well-defined functional on formal power 
series, valued in formal power series. Furthermore, all functionals ^ k and all operators 
Oj k appearing on the right side are taken at a reference point x. We now define the 
operator product coefficients by formula (|4*H|) in terms of the functionals ^f i ( . )j. Writing 
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out all terms explicitly, the interacting OPE-coefficients are thus given by 



:miH hm fc 

C n ...J( Xu Xn)j = ^(-l) fc — i 1 E 

' mi! -••mi.! 

fc=0 m,>l 1 re fe]<A 

^'(R mi (C», 1 (a; n );/^))^ {R m2 (O l2 (x n );I^))---¥^K mk (O lk+1 (x n ); J*™*)) 

/ n -n r \ 

x E II — ^(^r^t); I® 1 *) , (49) 



ni >0 \r=l 



where all local functionals ^ refer to the point x n . In order to make this formula well- 
defined, it is necessary to assume that the support of the cutoff function \ implicit in I is 
small enough so that the standard functionals are defined on the corresponding retarded 
product. However, this is no real restriction, because the OPE is an asymptotic short 
distance expansion, and we will later show that the coefficients do not depend on the 
particular choice of x asymptotically 

We claim that the coefficients Cj satisfy an OPE: 

Theorem 1. Let the interaction I = J Cdfibe renormalizable, i.e., [C] < 4. For a given 
tree T, 5 > 0, and given ii, . . . , i n G X, let 

A = 5+ (it[ij]J -depth(T), (50) 

and define the OPE coefficients C h „ An k by eq. (gnj). Then the OPE holds: 

O il (x 1 )iO i2 (x 2 )i ■ ■ ■ O in {x n )i ~ T ,<5 ^ C il i 2 ...i n k (x 1 ,x 2 , . . . ,x n )iOk{xn)i ■ (51) 

[fc]<A 



Remarks: 1) The theorem is false for non-renormalizable interactions. 
2) Since the topology on V ® J~(M, g) of which the interacting fields are elements is 
generated by a set of seminorms associated with functionals including the Hadamard 
states, it follows that the OPE will continue to hold on the factor algebra obtained by 
dividing by the Klein-Gordon equation, in the sense of expectation values in Hadamard 
states, to arbitrary orders in perturbation theory. 

Proof: Let Mi be defined as the remainder in the OPE, i.e., the left side of (|51jl minus the 
right side. We need to prove that e~ 8 Q(Nioip r (e)) tends to in the sense of distributions 
as e — > 0. The analysis of this limit is easiest in the case when T is the tree T = 
{S , Si, ... , S n } with one root S Q = {1, . . . , n} and n leaves Si = {i}. Then depth(T) = 
1, and ipriz) is the map that scales the Riemannian normal coordinates of the points 
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Xi, . . . , x n _i G U relative to y = x n by e, where U is a convex normal neighborhood of y. 
Thus, taking e = 2~ N , we must show that 

2 SN <S>(Af I (2~ N x 1 ,...,2- N x n )) asJV->oo, (52) 

in the sense of distributions valued in V, i.e., to any order in perturbation theory. In 
the above expression, and in the remainder of this proof, points Xj have been identified 
with their Riemannian normal coordinates around y = x n . In order to analyze the above 
expression, it is necessary to perform several intermediate decompositions of A/}, and we 
now explain how this is done. 

We first decompose Mi into contributions from the different orders in perturbation 
theory. The ring V of formal power series in the couplings contained in the interaction 
I = f £dfi can be decomposed into a direct sum 

V = ^^V(k), V(k) = Eigenspace of Kid/dKi for eigenvalue k, (53) 

k 

where the k-th summand corresponds to the k-th order in perturbation theory. Accord- 
ingly, A/} may be decomposed as ^2M( k ) into contributions from the various orders in 
perturbation theory, and likewise Cj = Yl^(k) e ^ c - Using the Bogoliubov formula, the 
k-th order perturbative contribution to A/} can be written in the form 

M (k) (x u ...,x n ) = J2 kl \' k .k n \ Tl Rk i {0{Xj)jmj) 

fciH Yk n =k j 

k 

- J2C {p )(x 1 ,...,x n )R k ^ p (0(x n ),I^ k -^), (54) 

p=0 

where the labels on O indicating the field species have been omitted to lighten the nota- 
tion. Let be defined as Af(k), but with the k-th order OPE coefficient C^) omitted. 
Then, using the definition of the OPE-coefficients, it can be seen that 

C il ... in J (x 1 ,...,x n ) {k) = W(A il ... in (x 1 ,...,x n ) {k) ) (55) 

and that M(k) = A(k) — Yl\j]<A ^K^k^Oj, where ^ are the free field reference functionals 
at point x n and where Oj are the free field Wick powers taken at point x n . Thus, the 
expectation value of the scaled, k-th order perturbative contribution to the remainder is 
given by 

$(M {k) (2- N x 1 , 2~ N x n )) = ^(A {k) (2- N x u 2- N x n )) 

- * J '(A w ( 2 " JV ^>--- ) 2 " JV ^))$(O j (2- Ar x n )), (56) 
b']<A 
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The right side of this equation is schematically of the form $(F) — $?y(F)$>(Ok(y)) , and 
for such expressions we will now write down an expression which will be useful to analyze 
the limit iV — > oo of eq. (|56jl . To derive this expression, perform a Taylor expansion with 
remainder about (y, . . . ,y) G U m of the m-th locally normal ordered product (|2()j) . 



m m 

■■Il*®>--- E n— 

*=1 |aiH h|o m |<p m «=1 

1 /-l m 

= - (1-tydp 1 :Jlmi)-n dt. (57) 
P- Jo i=1 

Here, the & denote Riemannian normal coordinates around y and are identified with 
points in R 4 , G N 4 is a multiindex, and quantities like |a$| are defined using standard 
multiindex conventions. As explained above, any element F G J~{M, g) supported in U 
may be written as a linear combination of expressions which consist of distributions u m 
supported in U m satisfying the wave front condition (jUJ), integrated with locally normal 
ordered products : 0® m : H . If we apply a Hadamard state $ to such an expression F, 
use the above Taylor series with remainder, and use the definition (|45)l for the standard 
functionals, then we get the following equation: 

$(F)- Y, * k y {F)®{O k {y)) 

[fc]<A 

= E E (A-m)! U ^---^^T'--C m 

m | a |=A-m+l V ' M 

pi m 

/ (l-t) A - m $(:a ftl #6)--^ m (te m ) :„)^A^). (58) 

i=l 

The key point to note about this identity is that there are now factors of £f 4 on the 
right side, which will work in our advantage when the points & are scaled by a small 
factor. On the other hand, the normal ordered expectation values in the second line are 
smooth (here we are using the assumption that $ is Hadamard), and so will not cause any 
trouble for such a scaling. We will now prove that (p)2*j) holds in the sense of distributions 
by exploiting this identity for F = An~\ in eq. (|56p. However, before we efficiently make 
use of that identity in (jHHj) . it is first necessary to rewrite A(k) in a suitable way, and to 
apply an induction in k. 

For this, we recall that the interaction Lagrangian density £ is confined to the convex 
normal neighborhood U since we are taking the couplings to be Ki(x) = Kix{x) with x a 
smooth cutoff function that is supported in U. We now "slice up" the support of C into 
contributions from different "shells" in U that are centered around y = x n , and that have 
thickness 2~ J , where j = 1, . . . , N. For this, we choose a compactly supported function •& 
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that is 1 on U, and we set 
Then C may be decomposed as 



&j{x) =<&(2 s x) 



(59) 



N 

C = -& N C + ^2('& j -i--& j )C. (60) 

3=1 

Each term in the sum is supported in a slice of thickness 2 _jf , see the following figure on 
p. 24. The key step is now to rewrite an interacting field quantity in a way that reflects 
the subdivision of the interaction region U into these slices. For this, we note that if Vj 
is the unitary in (|31|1 relating the interacting field with interactions Ij = J djCdp and 
Ij-i = f ■dj-xCdn, we have 

O^x) = V1V2 ■ ■ ■ VnO^^x)^ ■ ■ ■ V N )~ l , (61) 

for all x G U. Explicitly, Vj is given in terms of the relative S-matrix [3], 

k 

Here, p is any smooth function of compact support in U with the property that p(x) = 
for all x G J + (supp($i)) and p(x) = $o( x ) ~ fc> r all x G J~(supp(-#i)), and 

Pj(x) = p{2^x). Each term in this expansion can in turn be written in terms of retarded 
products Substituting the equation (|6Tj) into the formula for the remainder, and 
expanding in a perturbation expansion, we get the following identity: 

A (k) (2~ N x u 2-"x Tl ), = A {k) (2~ N x u 2- N x n ),„ 

+ E E II Sj Uto£) M 

P=0 k-p=k 1 -\ h k r 0<ati—<a r <N 

+ll + --- + lq 

■Af (p) (2- N x u ...,2- N x n ) lN 

■ n s s* P ASp* c )\k) ■ ^ 

0>Pl->Pq>N 

This complicated identity has the following structure. The sum on the right side is by 
definition only for p such that p < k, meaning that the terms under the sum only contain 
the remainder in the OPE up to (k — l)-th order in perturbation theory. This will enable us 
to use an inductive procedure to estimate the k-th order perturbative contribution to the 
remainder by the lower order contribution. The first term on the right side, A^){- ■ ■ )i N , 
is identical in nature with A^(. . . )/, with only exception that all the retarded products 
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implicit in its definition are now computed with respect to the interaction 1^ = j fi^Cdf! 
which is supported only in a small ball of radius 2~ N around y = x n . This will enable us 
to use a scaling argument to estimate this term. 

We now explain more precisely how the decomposition of Aha given in (|6l3|) will make it 
possible to analyze the scaling behavior of the fc-th order remainder in the OPE. For this, 
we take (Jo^j) . and substitute it into eq. (Joljj) . This gives an expression for $(A r (fc)(. ••)/), 
and to each term in this expression, we can apply eq. (p)%|) . Consider first the term arising 
from the first term on the right side of eq. ()63)1 . That term makes a contribution to 
$(jV(fc)(. . . )/) of the form 

$(A (fc) (2-V, . . . , 2~ N x n ) lN ) - * j {A {k) (2- N x 1 , . . . , 2- N x n ) lN m0 3 {2- N x n ) lN ) . 

b1<A 

(64) 

We must now substitute the expression for Aqa{. . . )i N . For this, we use that, on account 
of the Bogoliubov formula, the scaled interacting field with interaction 1^ is given by 
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0(2- N x) lN =^_^_ jf^R fc \ O(2- N x),0C(2- N yi )j HMVk)dfi N (y k ) , (65) 

where we have performed the change of integration variables yi — > 2 N y^ and where d[i^ is 
2 4iV -times the pull-back of d[i by the inverse of this map, which is smooth as N —>■ oo. If 
we now also use the Wick expansion of the retarded products (r3), along with the scaling 
degree property (r4) and combine the result with eq. (jSHJ), then we obtain that eq. 
scales as 2~~ N ( 5+1 \ as desired. 

Now we must take the second term on the right side in (JBlIj). substitute it into (foT)j) . 
and then analyze its scaling using (|58j). To do this, we must now proceed iteratively, in 
the order in perturbation theory k. For k = 0, there is nothing to show. For k > 0, we 
then inductively know the scaling (|52|) up to order k — 1, which enables us to estimate the 
remainders J\f(p\,p < k in the terms in the sum on the right side of (p)3"j) . More precisely, 
we may inductively assume that the p-th order remainder (p < k) has the structure 



^00(2-"*!, • • • , 2- N x n ) lN = 2-^+ A - +1 ) 



(A - to)! 

a|=A-m+l 



n m (2- N yi , . . . , 2- N y m , 2~ N x u . . . , 2"^)^ • ■ • 
/ (l-t) A ~ m : 8^(12-" yi )---d am (t2- N y m ): n dt /\dfi N ( yi ) , (66) 

where n m are the coefficients in a Wick-expansion of N( p )(xi, . . . , x n )i (note that we have 
also performed a change of integration variables yi — > 2 N yi). Using the scaling (r3) and 
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the fact that all terms A/"( p )(. . . )j may be written in terms of retarded products by means 
of the Bogoliubov formula, one can see that 

2 -^(b- 1 ]+-+W+m[£]-m) nm(2 -7V yi) ; 2- N y m , 2~ N xi, . . . , 2~ N x n ) ^ as N -> oo. (67) 

One now has to take formula l|66J). and substitute it into the sum on the right side of . 
From the product of Af( p )(- ■ -)i N with the relative ^-matrices there arise terms which 
blow up as N — > oo, and so these terms have to be carefully controlled. To understand in 
detail what type of diverging terms can arise, we must write out the explicit formula for 
the relative S'-matrices in terms of retarded products. Then we must write each retarded 
product in a Wick expansion (r4), and perform the products using Wick's theorem 
with uj replaced by H . Then we get a collection of terms, each of which is a product of 
H, r, n m and a locally normal ordered Wick power. These terms are evaluated on a set 
of spacetime arguments which are scaled by 2~ N ,2~ at , or 2~^\ and which are integrated 
against the compactly supported smooth functions $ or p. The arguments scaled by 
2~ a J arise from points in the interaction domain U within the ay-th slice, the arguments 
scaled by arise from points in the interaction domain U within the (3j-th slice, and 
the arguments scaled by 2^ correspond to the scaled arguments 2~ N Xi. More precisely, 
when we use Wick's theorem to perform the products in the second term in (|63|h there 
arise "contractions" between points in the a«-th and /3^-th slice, indicated by lines in the 
following figure: 




Each such contraction is associated with a factor H(2~ a iyi,2~ l3i y 2 ) (or a derivative 
thereof), which, using the explicit form of the Hadamard parametrix H, is seen to scale 
as 2 2mm ( Qi ' Q j), in the sense of the scaling degree of a distribution (with a correspondingly 
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larger power when derivatives are present). Furthermore, the scaling of the retarded 
products in a term in eq. associated with the i-th slice may also be controlled. 
Namely, using the Wick expansion (r3), we see that a retarded product associated with 
the i-th slice contributes factors of r(2~ l yi, . . . ,2~ l yi), the scaling power of which may 
then be controlled using (r4). Finally, the scaling of n m (2~ N xi, . . . ,2~ N y m ) is controlled 
by eq. (|67|). Thus, the rate at which the terms in the sum on the right side of (|63|) blow 
up can be controlled. We finally need to substitute each such term into eq. (|5l)|). and 
use (I58J) . If we carefully keep track of all the scaling powers, then we find that a typical 
term contributing to N(k) arising from these substitutions has the scaling power 

2 N{-5-Ak-l)+[C] £V kjaj+lC] Ei hPi /gg) 

Using a geometric series, the sum of such terms is estimated by 
fc-i 



P=0 k — p= fei H h k r 1 < m ■ ■ ■ < a r < N 

-Hi H \-lg K /3i • • • < Pq < N 



< const. 2 - N{ - s+l -^-^ < const. 2 -^ 5+1 \ (69) 



where we have used in the last step that the interaction is renormalizable, [£] < 4. Thus, 
the total scaling of the sum of terms in J\f^)(2~ N xi, . . . , 2~ N x n ) is given by 2~ Ar ( <5+1 ), which 
implies the convergence of fl52(l . On the other hand, for non- renormalizable interactions, 
we would not get convergence. 

The analysis for a general tree T is in principle not very different from the one just 
given. For a general tree, the minimum distance between points in a scaled configuration 
x{2~ N ) is of order 2~ depth ^ N , and not just 2~ N as in the simple tree studied above. This 
implies that the scaling of the corresponding quantities in the perturbative expansions is 
different. One now has to go through the above steps again and take that different scaling 
into account. If this is done, then the result claimed in the theorem is obtained. □ 



4 Properties of the OPE coefficients 

We would now like to establish a number of important general properties of the OPE 
coefficients defined in the previous section. These properties are 

1. Microlocal spectrum condition. 

2. Locality and covariance. 

3. Renormalization group. 

4. Associativity. 
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5. Scaling expansion. 



Except for the last one, these properties were postulated as axioms in the paper so 
our present work can be viewed as a confirmation of |25j . 

We first establish the microlocal spectrum condition. With the graph theoretical 
notation Q n<m introduced in section 2, let us define the following subset of the cotangent 
space T*ikP: 



F n ,m(M,g) = jOi, h; . . . ; x n , k n ) G T*M n \ {0}; 3 decorated graph G(x,y,p) G G m , n 
such that hi = p e — p e for all Xi and 

e:s(e)=Xi e:t(e)=Xi 

such that = p e — p e for all yi and 

e:s{e)=yi e:t(e)=yi 

such that yi G J + ({xi, . . . , x n }) fl J~({xi, . . . , for all 1 < i < m, \ . 

(70) 

The microlocal spectrum condition for the OPE-coefficients is statement that 



WF(C) \U n C |J r -,n. (71) 



m>0 



where U n is some neighborhood of {{x,x, . . . ,x) G M n }, and where "WF" is the wave 
front set of a distribution |2*7j . 

The microlocal condition in the above form (|7T|) is similar in nature to a condition that 
was obtained by [3] for the n-point correlation functions of Wick powers in Hadamard 
states in the context of linear field theory. The difference to the above condition is that also 
interaction vertices are now allowed, which were not considered in jSj. These interaction 
vertices correspond to the contributions m > 1 in (|71|) and genuinely weaken the bound 
on the wave front set relative to the linear case (for n > 4). The interaction vertices 
arise from the non-linear interactions present in the theory. As we will see, the maximum 
valence of the interaction vertices allowed in WF(C) is equal to the maximum power of 
the field <fi that appears in the interaction Lagrangian C. Since we restrict ourselves to 
renormalizable interactions in 4 spacetime dimensions, that maximum valence is equal to 
4. Note, however, that the wave front condition (|71jl is only an upper bound, and does 
not say whether interaction vertices will actually contribute to WF(C) or not. We have 
checked this for the OPE-coefficient in front of the identity operator in the expansion 
of the product 4>(xi)i ■ • • <f>(x^)j of four interacting fields, to first order in perturbation 
theory in Minkowski space, where a contribution from an interaction vertex in T^i would 
be allowed according to the above estimate (fTTf . Using our definition pHj) of C and using 
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the integrals in Pj, we found in this example that such a contribution is actually absent 
from WF(C). Hence, the estimate (j71|) is not sharp. 

Let us now prove the microlocal condition ([7T| . By eq. (|49j) . the microlocal spectrum 
condition will follow if we can show that if 



then WF(m„) C U m r n m . To prove this statement, we expand the retarded products as 
in eq. (|27j). then multiply them using the Wick expansion formula pOjl . and finally apply 
the functional ^ . The result will be sum of expressions each of which is a product of 
r's, of H 's and expectation values in of locally normal ordered expressions, which are 
integrated over interaction vertices against the smooth test function \ of compact support 
appearing as infrared cutoff in the Lagrange density, C. The expectation value in ^ of any 
locally normal ordered expression is smooth, the wave front set of the r's is given above 
in (r6), see (jSHJ, while the wave front property of H is WF(H) = {(yi, k±; y^ fe); k\ = 
p, k 2 = —p, p G dV+}, where p is a coparallel, cotangent vector field along a null geodesic 
(edge) joining 2/1,2/2- We now combine these facts using the wave- front set calculus of 
Hormander, by which we mean the following theorems about the behavior of the wave 
front set under the operations of smoothing, and products [2Zj : Let X, Y be manifolds 
(in our applications, they are Cartesian powers of M). If K G D'(X x Y) is a distribution 
and / G T>(Y) a smooth test function, then the distribution u(x) = f Y K(x, y)f(y)d/J l (y) 
has wave front set 



Secondly, let u, v G V(X) so that [WF(u) + WF(v)] n {0} = 0. Then the distributional 
product uv is defined and has wave front set 

WF(m>) c {(x,k + p) G T*X; (x, k) G WF(«) U {0}, (x,p) G WF(t>) U {0}} . (74) 

Applying these rules to the above products of r's and if 's, we essentially obtain that 
WF(u n ) is a subset of U m r nim . For example, the momentum conservation rule in the 
third line of eq. (fTHjl follows from the additive and smoothing properties (fTHjl . (|74*j) com- 
bined with the fact that we are integrating the interaction vertices against the smooth 
testfunction, x, of compact support. Similarly, we obtain the second line from the additive 
property ([74)1 . Finally, we need to prove the support restriction on the interaction vertices 
in the fourth line of (|70|) . For this, we note that the contribution to the wave front set 
of u n from the interaction vertices y^ arises only from points that are in the support of 
the interaction, 7, i.e., in the support of %■ U be an arbitrary small neighborhood 
of J~({xi, . . . , x n }) PI J + ({xi, . . . , £„}), and let x' be a cutoff function supported in U. 




(72) 



WF (u) C {{x, k) G T*X; (x, k; y, 0) G WF(K)} 



(73) 
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Then x — x' is supported outside of the domain of dependence D({xi, . . . ,x n }), and so 
there exists by (pTTj) a unitary V such that 

J] O lk (x h ) P = V I J] O^x.h \v* . (75) 

k I k ) 

Thus, because of (|IHJ), we see that changing x to x is equivalent to changing the standard 
functionals from . . ) to ^'(V . . . V*). We claim that this would not, however, change 
our above wave front argument. Indeed, the only property of the functionals that was used 
in the above wave front set argument was that the expectation values of locally normal 
ordered expressions in ^ are smooth. This does not change if we change the standard 
functionals from . . ) to ^ 3 (V . . . V*). Consequently, we have shown that contributions 
to ()7()jl arise only from interaction vertices y& in U. Since U was an arbitrarily small 
neighborhood of J~({xi, . . . , x n }) fl J + ({xi, . . . , x n }), the support restriction in the last 
line of (|2J|> follows. 

We next show that the OPE coefficients have the following local and covariance prop- 
erty: Let / : M — > M' be a causality preserving isometric embedding, let Cj respectively 
Cj, be the OPE coefficients on the respective spacetimes, and let 5 > be given. Finally, 
assume that there are open neighborhoods U C M and U' C M' with f(U) C U' where 
the cutoff functions x respectively x' implicit in the interactions / and I' are equal to 1. 
Then, supposing that A in eq. ({H| is chosen as in Theorem 1, we have on U 

f*C' r ~ Ti5 Cj for all trees T. (76) 

This condition essentially follows from the fact that the interacting fields are local and 
covariant, in the sense of (jlHj) . which follows in turn from the fact that the individual terms 
in the perturbation expansion of the interacting fields are local and covariant. However, a 
complication arises from the fact that the algebra embedding af in (JT3j) is not simply given 
in terms of the corresponding free field homomorphism, but is more complicated • 

Instead of taking into account the more complicated definition of aj at the interacting 
level, one can also more directly prove (|76|). For this, we note that, if the cutoff function x' 
on M' were such that f*\' = X; then we would have equality in (fTHjl. because the retarded 
products and standard functionals which are the ingredients in the definition of Cj have 
a local and covariant dependence simultaneously on both \ and the metric g implicit in 
J, by property (r5). Thus, it is sufficient to show that Cj is essentially independent of 
the cutoff function x- in other words, if x and x' are two cutoff functions (on the same 
spacetime) which are equal to 1 on U, and if Cj and Cr are the corresponding OPE 
coefficients, then we must show that 

C/ ~<5,r Cp (77) 

holds on U. To prove this statement, we simply apply the functionals *ff(V . . . V*)j to the 
remainder Mp of the OPE formed using the coefficients Cp, where V is the unitary in 
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eq. (|3Tj) relating the interactions I and V . Then we find ty(VJ\fpV*)i = Cj — Cp. Since 
A/j' is the remainder of the OPE, and since ty(V . . . V*)i is a Hadamard functional, it 
follows that e~ s (Ci — Cj/) o ^rie) will go to zero by theorem 1, which is what we needed 
to show. 

In [23! 5 ft was shown that the perturbative interacting fields obey a "local covariant 
renormalization group flow". The construction of this flow involves the consideration of 
a 1-parameter family A 2 g of conformally rescaled metrics, where A £ H>cb and states 
how the interacting fields Oi(x)i change under such a rescaling. In [23], a simple general 
argument was given that the existence of such a local covariant renormalization group 
flow implies a corresponding flow of the OPE coefficients if the theory possesses an OPE 
in the sense described in the previous section. The key assumption on the nature of the 
RG made in [23] was that there exists a suitable "basis" of functionals which is in some 
sense "dual" to the fields. This is the case in perturbation theory, on account of (JEJ). 
Hence, it follows by the argument of [25j that 

Z(\y- n ■■■Z(\) i \[ t Z(\)- 1 ] k l C il ... in k [M,g]j C n ..,J[M,X 2 g} Iix) (78) 

for all trees T. Here, the prefactors are linear maps (whose construction and properties 
was described in [23]) 

Z(\) i j £ End(Ej, Ej) , (79) 

where Ej is the vector bundle in which the field Oi lives, and /(A) is an interaction of 
the same form as / with suitable "running" couplings ftj(A), whose construction was also 
described in |23j . 

An associativity property for the OPE coefficients may be formulated as follows 
(see j25J for details). Let x(e) be a curve in configuration space representing the merger 
of the points according to a tree T, that is, 4>r(£) '■ x h- > x(e). In this situation, we should 
be able to perform the OPE successively, in the hierarchical order represented by the tree, 
thus leading to some kind of "asymptotic factorization" . That is, we should be allowed to 
first perform the OPE for each subtree, and then successively the OPE's corresponding 
to the branches relating the subtrees, and so fourth. For example, for the tree T given in 
the figure on p. 16, we should be allowed to perform the OPE successively as indicated 
by the brackets {Oi02){0^0 ^) . To formulate this condition more precisely, we recall the 
notation s(e),t(e) for the source and targets of an edge, e, in the tree T. Furthermore, 
for x = (x±, . . . , x n ) £ M n , and for each node S £ T of the tree, let us set 

xs = x m (s), m(S) = max{i : i £ S} . (80) 

Finally, we consider maps i : T — ► X which associates with every node S £ T of the tree 
an element is £ X, the index set labelling the fields. With these notations in place, the 
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associativity property can be stated as follows. Let 5 > 0, let T be a tree, and let Mg be 
the set of all "spacelike configurations" x — {x\, . . . , x n ) 

= {xe M"; Xi <£ J + (xj) U J~(xj) forall i, j} . (81) 

Then, on Mg, we have 

C^^.i^ (yXl, . . . , Xn)j 

~ T > S II C i*tW> e such that s ( £ ) = 5 >' S ({^(e)! e SUch tllat S ( e ) = ' ( 82 ) 

where the sum is over all z, with the properties that 

i{k} = k, k = l,...,n, i{i,..., n } = j, (83) 
[it(e))<S s VSeT, (84) 

e:s(e)=S 

where 5s > are chosen sufficiently large. Note that it makes sense to consider the relation 
~s,t with respect to the open subset of spacelike configurations, because a configuration 
remains spacelike when scaled down by ipr(e), at least provided the points x G Mg are in 
a sufficiently small neighborhood of the total diagonal, which we assume is the case. The 
technical reason for restricting the OPE to pairwise spacelike related points is that the 
OPE coefficients C are smooth on Mg, by the microlocal spectrum condition (see eq. ()71jl ). 
and so convergence in the sense of ~s,T is more straightforward to study. Furthermore, 
since the interacting fields commute for spacelike related points, there are no ordering 
issues when working with the configurations in Mg. Also, from a physical viewpoint, the 
notion of "short distances" is somewhat unclear if lightlike directions are included. 

In it is shown that associativity in the above sense is an automatic consequence 
if the theory also possess a suitable local covariant renormalization group with suitable 
properties, and if the OPE holds not only for each fixed spacetime and fixed choice of 
couplings, but instead also uniformly in a suitable sense for smooth families of metrics and 
couplings (termed "condition (L)" in that paper). As we have already described, the ex- 
istence of a local covariant renormalization group in perturbation theory was established 
in [2B]- It is a general property of the perturbative renormalization group that Z l j(X) is 
given by A~W times a polynomial in In A at each finite order in perturbation theory j2H|, 
and that the running couplings K{(\) in J(A) have a power law dependence A 4- ^, which 
is modified by polynomials in In A at any given order in perturbation theory. These are 
essentially the properties for the renormalization group required in (L) at the pertur- 
bative level, except that the running of couplings /(A) is not exactly smooth at A = 
as required in (L), but instead contain logarithmic terms at each order in perturbation 
theory. However, the argument given in [22] is insensitive to such logarithmic corrections. 
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To also establish the desired smooth dependence of the OPE-coefficients under smooth 
variations of the metric required in (L) at the perturbative level, it is necessary to go 
through the proof of Theorem 1 for families of spacetimes and corresponding families of 
states depending smoothly on a parameter in the sense [22] arid analyze the behavior 
of the constructions under variations of the parameter. This can indeed be done, using 
the smooth dependence of the retarded products under such parameters [22], as well as 
the techniques and type of arguments employed in the appendix of [24 . However, even 
though the repetition of these arguments is in principle straightforward, that analysis is 
quite lengthy and cumbersome, and not very illuminating. It is therefore omitted. Since 
a perturbative version of condition (L) holds, Theorem 1 of [2H] then implies that the 
associativity property holds on the space Mg of pairwise spacelike configurations. 

In the remainder of this section we will prove that the OPE coefficients themselves can 
be expanded for asymptotically small distances in terms of curvature terms and Minkowski 
distributions in the tangent space (for spacelike related configurations in Mg, to which 
we shall restrict ourselves in the remainder of this section). The construction of this 
"scaling expansion" involves the Mellin transform, Ai[f,z], of a function f(x) defined on 
IR>o vanishing near infinity, with at most polynomial type singularity j3H] as x — » 0. It is 
defined by 

POO 

M[f,z]= / x iz - 1 f{x)dx, (85) 
Jo 

and is an analytic function of z for sufficiently small lm(z) < y where y depends upon 
the strength of the singularity of /. The inverse Mellin transform of F(z) = M[f, z] is 
given by 

I r+ioo+c 

M~ 1 [F,x] = — x z F(z)dz (86) 

2 71 "! J -ioo+c 

where the integration contour is to the right of all poles of F(z) in the complex z-plane. 
The Mellin transform is useful in the context of functions f(x) possessing near x = an 
asymptotic expansion of the form 

a Pil x-P(\nx) 1 , (87) 

p l 

where p is bounded from above, and where the sum over I is finite for any p. It can 
be seen that the Mellin transform of such a function possesses isolated poles at z = ip 
in the complex plane, with finite multiplicities. Furthermore, the asymptotic expansion 
coefficients a Pt i are the residues of the Mellin transform, i.e., 

a P ,i = ^Res 2= _ ip |(^ + ip) l M[f, z}\. (88) 

We now define, for each tree T, distributions Cj that give the desired scaling expansion 
of the OPE coefficient Cj relative to the scaling function ipri^) '■ Mq —>■ Mq, x t— > x(s) 
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defined above in eq. (j39j). by extracting the poles of Cj o ipr{^) in £ using the Mellin 
transform. In order to do this, let x = (xi, . . . ,x n ) G Mg be a spacelike configuration 
of n points, let T be a tree, and let e i— > be the corresponding curve in Ms*. By 
the microlocal spectrum property, the OPE coefficients are smooth on Mg, so we may 
consider [Cj o i/)t(£)](x) = Ci(x(e)) as a smooth function in e at any fixed value of the 
argument x. As we will show in the proof of the next theorem, if we fix the parameter 
5 > in the operator product expansion, this function has an expansion of the form 

C 7 (f(e)) =££M*>-(ln £ y + . . . , (89) 
p * 

near e = 0, where the dots stand for a remainder vanishing faster than e s . Here, p is in 
the range from —5 to D = depth(T) ■ (—[k] + ^2j[ij])-> an d the sum over I is finite for each 
p, at any given order in perturbation theory 7 . Consequently, we can define the Mellin 
transform 8 of this function in the variable e 

POO 

M r (x, z) =M\C!o ^ T (e), z] = Ci{x{e)) e lz ~ l de , (90) 

Jo 

which is now a function of x £ that is analytic in z G C for sufficiently small lm(z). 
Furthermore, by the above expansion ([89)1 . it is meromorphic on a domain including 
lm(z) < 6, with poles possibly at iS, i(6 — 1), i(8 — 2), . . . , —iD. Let us now choose a 
contour C around these points as illustrated in the figure. 
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Define 

C r (f) 7 = — / M T (x , z) dz . (91) 
2vn J c 

Concerning this function on Mg, we have the following theorem. 
7 Note however that the range of I increases with the perturbation order. 

To make this expression well defined, we need to arbitrarily cut off the integral for large e (where the 
map ipr(£) is not well-defined anyway). How this cutoff is chosen does not affect the following discussion. 
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Theorem 2. 1. We have Cf ~r,5 Ci for spacelike configurations, and therefore 
O il (x 1 ) I Oi 2 (x 2 )i ■ ■ ■ O in {x n )i ~ r ,5 

[fe]<A 

on Mg. 

2. Cj is local and covariant, i.e., if / : (M,g) — > (M',g') is an orientation, causality 
preserving isometric embedding then Cj = f*Cjf . In particular, Cj does not 
depend on the choice of cutoff function \ used in the definition of the interacting 
field. 

3. The expression Cj{x) is the sum of residue of M r (z, x) corresponding to the poles 
in the contour C in (JHTJ), 

Cj(x) = £ Res 2= _ ip {M T (x, z)} . (93) 

p>-6 

These have the following form: 
Res z= . ip {M T (x, z)} = J2 p aW iai ■ ■ ■ V Qfe) iW 2 „3M4 0n)] W a (^, . . . , Cn-l) , (94) 

a 

where P a is a polynomial in the Riemann tensor and (finitely many) of its covari- 
ant derivatives evaluated at x n , valued in some tensor power of the tangent space 
(T Xn M)® a , while & are the Riemannian normal coordinates of X\, . . . ,x n _i around 
x n , identified with vectors in K 4 via a tetrad. The sum over a is finite, and each W a 
is a Lorentz covariant distribution on IR 4(n-1 ) (defined on spacelike configurations), 
that is valued in (R 4 )® a (identified with (T* n M)® a via the tetrad), depending poly- 
nomially on m 2 , a. Thus, for any proper, orthochronous Lorentz transformation A, 
we have 

W a (A^ . . . , A£n-i) = D b a (A) W 6 (6, . . . , Cn-i) , (95) 
where -D(A) is the corresponding tensor representation. 

4. There are Lorentz invariant distributions V a i such that 



N 



(96) 



where the sum is only over a finite range of /, at any given, but finite order in 
perturbation theory. 



33 



Remarks: 1) Equations ([93)1 . (}9*4*|) constitute the claimed scaling expansion. That scal- 
ing expansion is similar in nature to a corresponding expansion derived in [22] for the 
short distance behavior of time ordered products. However, note that the scaling expan- 
sion above is more general than that derived in [22] , because it involves the consideration 
of more general scaling functions ipr{^) corresponding to general trees. Also, the Mellin 
transformation technique was not used in [22J. 

2) The restriction to spacelike configurations has mainly been made for technical 
reasons, to avoid technical issues that could arise when taking the Mellin transformation 
of a distribution. However, we expect that all constructions and properties summarized 
in the above theorem hold for all configurations, i.e., in the sense of distributions on M n . 

Proof. Let us first argue that the claimed meromorphicity of the Mellin transform 
of C(x(e))i holds, or equivalently, that it has an asymptotic expansion as claimed in 
eq. (JHH)- For this, we first consider the simple tree T = {S , Si, ... , S n } with one root 
So — {1, • • • , n} and n leaves Si = {i}. The depth of this tree is depth(T) = 1, and ipri^) 
is simply the map which multiplies the Riemann normal coordinates of Xi relative to x n 
by e. To analyze the corresponding scaling of the OPE coefficients, we note that, by the 
local and covariance property of the OPE coefficients, a rescaling of the arguments Xj is 
equivalent to changing the metric from g to s* £ g, where s £ : M — ► M is the diffeomorphism 
that scales the Riemannian coordinates of a point around x n by e. Thus, we have have 

Cj[M, g] o i>? (e) ~ T)(5 d[M, s* £ g] (97) 

for this tree. Next, we use the fact that in perturbation theory, there exists a local and 
covariant renormalization group [23! • This implies that, up to terms of order e s 

Z(e)\ ■ ■ ■ Z(e)\ i'ZieyVC^iM, g(e)} m = C n ..,J[M, g]j o ^ r (e) , (98) 

where g(e) = s 2 s*g. The point is now that the metric g(e) has a smooth dependence upon 
e, as may be seen by rewriting it in Riemannian normal coordinates as g llv {e£ s )d^d£ s v . 
The OPE coefficients in turn have a smooth dependence upon smooth variations of the 
metric, since all the quantities in their definition have this property [22] • The only singular 
terms (in e) in the expression on the left side can therefore come from (a) the running 
couplings in 1(e) and (b) the Z(e)-factors. However, by the general analysis of the local 
renormalization given in [23], these quantities are polynomials in e~ p (lne) 1 of finite degree 
to any given order in perturbation theory. Thus, up to terms vanishing faster than e s , 
the quantity Cj(x(e)) has an expansion of the type (jH9*|l . as desired. 

The argument just given may be generalized to arbitrary trees by an induction in the 
depth of the tree. For trees of depth one we have just proven the statement. Let us 
inductively assume that we have proven ([H9*|) for trees of depth d. To deal with trees of 
depth d+1, one notices that a tree of depth > 2 can always be decomposed into a tree 
S of depth one connected to the root, and trees 7{, . . . ,% attached to the leaves of S. 
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Thus, we may write 

r 

T = SU{J% 

see the figure for an example with r = 2. 



(99) 



<S — {So, Si, S2} 



— \S\, S3, 5*4} 




T~l — {S2, S5, . . . , Sg} 



The key point is now that the map V't(^) factorizes by the inductive nature of its 
definition (}39|) . while the OPE coefficient Cj factorizes by the associativity property (J82|) 
under this decomposition of the tree. This gives 



C h ... in k o <4> r (e) 



J T,S 



SgLeaves 



(T t )} n °ip%{e) 



t=l 



(100) 



where we note that, since S has depth one, the map ips{ e ) is given in terms of the diffeo- 
morphism s £ : M — » M which scales the Riemann normal coordinates around x n of points 
by e. On the right side of this expression, we can now apply the induction hypothesis 
to the expression in brackets, because each of the trees % has depth < d. Furthermore, 
since ips{^) is given by the diffeomorphism s £ , we may again apply the general covariance 
and renormalization group property as we did above to convert the action of s £ into a 
smooth change g(e) in the metric and Z(e)-factors depending polynomially on e~ p (lne) 1 . 
This proves the equation (}89|) . 

We next come to the actual proof of the theorem. To prove 1), we need to show that 
e~ s (Cj — Cj) o ipri.^) — * as e — > 0, to any finite but arbitrary order in perturbation 
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theory. We have, using the definition of Cj and the relation ipri^) ° i'ri^') = tpri^') 

POO 

Cjo^ T (e) = ^Res 2= _ ip / (Cj o ^ r {e')) o ^( e ) e'^ 1 de' 

P>-6 J ° 



p>-8 ^ ' 



= ^e- p ^^Res, = _ ip {(z + ip)^ r (f,z)} 

p>-S i 

= EE £ ~ p ( ln£ )Vi(^)» (ioi) 

p>-6 I 

where we have performed a change of integration variables in the second step. Comparing 
with (|89p. this formula implies that C 7 (x(e))i differs from C(x(e))j by a term vanishing 
faster than e s . This proves the assertion 1). 

To prove 2), we recall that we have already proven above that, at each order in 
perturbation theory, f*Cp[M', g'\ ~s,t Cj[M, g], so the difference between these two terms 
vanishes faster than e s . This difference term will change the Mellin transform Ai T (z) only 
by a term that is analytic in a domain including lm(z) < 5, and thus will not contribute to 
Cj[M,g) respectively Cj,[M' , g'], because the contour integral of a holomorphic function 
vanishes. 

For 3), consider again the metric g(e) = e 2 s*g. Its components in Riemann normal co- 
ordinates around x n have a Taylor expansion of the form g^ u (et;) = £^Pv(V fc -R aia2a3a4 (x n ), 
where P/v are polynomials. Since the Cj are local and covariant by 2), it follows that 
they can be viewed as functionals of the Riemann tensor and its covariant derivatives at 
point x (or £ = 0) which enters via P^. Thus, it follows that the iV-th e-derivative is 



aN 

( Cj[g(e)](x)\ e=0 = wNl --- Nr (0\{ P NA^^---V ak) R m ^n),C 3 



8e N 

N 1 + ---+N r =N 



where 



r) r C r 



(102) 



(103) 



dP Nl ---dP Nr 

We define W a and P a by the above relation, i.e., P a is the appropriate product of the Pn v 
with the polynomial ^-dependence taken out and absorbed in the definition of W a . To 
prove the desired relation (|51jl. we must now show that the ^-derivatives of Cj[g(e)] at e — 
vanish for N sufficiently large. This quantity arises from the quantity Ci[g(e)} o ipr{^') 
by taking a Mellin transform in e', then taking N derivatives with respect to e at e = 0, 
and finally extracting the residue in z. It follows that we only need to show that the 
quantity obtained by taking N derivatives with respect to e of Ci[g(e)} o ipq-(e') vanishes 
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faster than e' 5 , because such a term would not give rise to poles of the Mellin transform 
in the domain Im(z) < 5. Consider first the case when T has depth one, so that ipri^') 
is simply a dilation of the Riemann normal coordinates by e' . As above in eq. ()98)) . by 
combining the renormalization group and general covariance, the action of such a dilation 
may be translated into changing g(e) to g{ee f ), along with a suitable set of Z(V)-factors, 
and running couplings in I(e'). The point is now that the N derivatives on e will produce, 
when acting on g(ee'), precisely N positive powers of s'. Thus, if iV is sufficiently large, 
then the resulting positive powers will dominate the corresponding negative powers in 
the Z(£')-factors, and we get the desired result. The generalization of this argument 
to arbitrary trees T can be done via an induction argument based upon formula 
similar to the one given there. 

To prove the Lorentz invariance of the W a , let us now define, for each Lorentz trans- 
formation A G SO(3, 1)J the diffeomorphism f\ : £ i— > A£. It defines a causality and orien- 
tation preserving isometric embedding between the spacetimes (M, g) and (M, f^g = g\) 
with the same orientations. Thus, using the local covariance property f\C r [g] = C r [g\], 
and the transformation property of f^P" 1 = D a b(A)P b under this diffeomorphism, it fol- 
lows that 

J2PaW k R aia2a3a4 (x n ),tn W a (ACi, . . . , Af n _i) 

a 

= Pa[V k R aia2a3a4 (x n ),a] D a b (A)W b {^, . . . , £ n _0 . (104) 

a.b 

However, since this holds for all metrics, eq. follows. 

Using the definition of W a just given, item 4) immediately follows from the fact that 
Cj has an expansion of the form (JHHJ), and that Cj ~ sr C v □ 

5 Example 

We now illustrate our general method for computing the OPE in curved spacetime by an 
example. Let us summarize again the steps needed in this computation. 

1. Fix a desired accuracy, S, of the OPE, and determine A as in Theorem 1. 

2. Identify the retarded products in eq. (jH^ that are needed to compute the desired 
OPE coefficient to a given order in perturbation theory and a given accuracy 5, and 
determine them, using the methods of the papers fH\ |2"2"] . 

3. Perform a local "Wick expansion (|27j) of all retarded products. In places in (J49J) 
where two retarded products are multiplied, perform the products using Wick's 
theorem (flT)|) (with lo^ in that formula replaced by H). Apply the standard func- 
tional (jH]]) to the resulting expressions in the way indicated in (I49J) . 
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This yields the desired OPE coefficient. If one is interested in the asymptotic behavior 
of the OPE coefficient (up to order 5 in e) under a rescaling of a point ipr{^) '■ % x(e) 
associated with a given tree T, perform the following step: 

4. Take the Mellin-transform of Ci(x(e)) in e as in eq. (HIIJ), and define Cj(x) to be 
the sum of its residue at the poles id, i(8 — 1), . . . , as in eq. (jHIjl . The result then 
automatically has the form of curvature terms times Minkowski distributions in the 
relative coordinates as described in item 2) of Theorem 2, and it is equivalent to Cj 
in the sense that Cj ~r,<5 C l . 

As an illustration, we now apply this method to the determine the coefficient Cj in 
the triple product of operators 

0(xi)/0(x 2 )/0(x 3 )/ = h C(xx, x 2 , x 3 )i (f)(x 3 )i + ... (105) 

up to first order in perturbation theory in the interaction / = f M Cdfi — — 4 J M K(p A d[i, 
and accuracy 5 = 0. We then discuss the scaling expansion as in item 2) of Theorem 2. 
As discussed above, we impose an infrared cutoff by taking k(x) = nx( x ) & t intermediate 
steps, where x is a smooth cutoff function, but the final answers will not depend on the 
choice of x, see eq. (ITTj) . The different ways of scaling the 3 points together give rise 
to different limiting behavior Cj of Cj. We choose to investigate the most interesting 
case when all points are scaled together at the same rate, i.e., under the scaling map 
V't(s) : x — > x{e), where T = {{1,2,3}, {1}, {2}, {3}}. If are the Riemannian normal 
coordinates of the points %i around x 3 , this corresponds to £(e) = e&i ££3) with 
6 = 0. 

Consider first the zeroth order perturbative contribution. According to eq. ()49j1 . this 
is given for a general OPE coefficient of a triple product of operators by 

C il i 2i3 j (x 1 ,x 2 ,x 3 ) (0) = ^(C n (x 1 )C i2 (x 2 )C i3 (a;3)), (106) 

where the reference point for the functional ^ (see eq. ([45))) is x 3 throughout. We are 
interested in the case O ix = O i2 = O ia = Oj = <fi. In order to determine the action of the 
functional ^l/- 7 in this case, we need to perform the local Wick expansion of the product 
0(xi)0(x2)0(x3), which is given by 

4>(xi)4>(x 2 )(f)(x 3 ) = : 0® 3 : H (xt, x 2 , x 3 ) + H(x 1 ,x 2 )4>(x 3 ) + cyclic(l, 2, 3) (107) 

Applying now the definition of the functional ^ (with Oj chosen to be <ft, and reference 
point x 3 ) gives 

C(x!,x 2 , x 3 )( ) = H(xi, x 2 ) + H(x 2 , x 3 ) + H(x h x 3 ) , (108) 

for any S. 
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In order to determine the representer Cj to zeroth order in perturbation theory for 
our choice 5 = (see eq. (jM}), we are instructed to compose the above result with the 
map ipr{e), then take a Mellin transform in the variable e, and then extract the residue 
at the poles 0, — i, — 2i in the complex z-plane via the contour integral As explained 
above, taking a Mellin transform and then extracting those residue is a way to compute 
the corresponding coefficients of e°,£: _1 ,e _2 in the expansion ()89j) of the distribution 
Cj o ipr{s) in e. In the present simple example, it is easier to compute the coefficients 
directly from eq. (jl08j) . by using the corresponding short distance expansions [TH! around 
£3 of the quantities a, v n appearing in the Hadamard parametrix H. The result is 



^(x 1 ,x 2 ,x 3 ) (0) =^x:{(^^ 

-iHn(fc-&) a }, (109) 



3 (6 - e,) 4 " 6 (6 - £,) 2 12 



where all curvature tensors are taken at 23. By item 1) of Theorem 2, we know that 

Cj ~t,0 Cj". 

We next consider the first order perturbative contribution to a general OPE coefficient 
for a general triple product. By formula (|%9"j) this is given by 



C il j 2 i 3 : '(Xi,X2,X3)(i) = i 



M 



^(O n (x 1 )^ 2 (x 2 )R 1 (O l3 (x 3 ),/:(z/))) + cyclic(l,2,3) 

X(y)dfi(y). (110) 



- £ ^(R 1 (O fc (x 3 ),/:(y))vI/ fc (O il (x 1 )O i2 (x 2 )O l 3(a;3)) 

[fc]<A 

We are again interested in the case = Oi 2 = Oi 3 = Oj = 4>. The constant A depends 
on the desired precision of the OPE governed by 5, and is given in Thm. 1. For our choice 
5 = 0, we have we have to choose A = 3. In this case, it can be seen that only Ok = 3 
will make a contribution. Thus, the required retarded products and their Wick expansion 
in eq. ()110|) are 



Mn^riv)) = E ( n - j b"i^- jb )i r *^y) : ^' k ^)r-\y) * (m) 



for n = 1, 3 and m = 4. Consequently, we need to know for k = 1, 3. The method [22] 
for defining the "renormalized" distribution gives r x = iAjj, where Ajj is the retarded 



39 



propagator, and 

rs = -^A'fa-tfEifllM] (ii2) 



32 u V <T + iOt 

3 2 v- o n ( lnV + iOt)+jjn(a + iOt)' 

+ -Vn > V n +l 7 A 6 — t) - — 

n \ 

~, ^ v n+m ( ln 2 ((J + i0t) ^ a k , „ ,\ , 

m,n \ k / 

where t(x,y) = t(x) — r(y), for any time function r : M — > R. The symbol A stands for 
the operator 9 

A = D + (V At lnD)V M , (113) 

where D is the VanVleck determinant defined above in eq. (fTTj) . The desired OPE coeffi- 
cient can now be obtained. The result is 



C(x l ,x 2} x 3 ) {1) =in / [H(xi,y)H(x 2 ,y)r l (x 3 ,y) 



M 



+ cyclic(l, 2, 3) - r 3 (x 3 , y)]x{y) dfx(y) , (114) 



where H is the local Hadamard parametrix. To get the desired representer CJ provided 
by Theorem 2, we again use the definitions eqs. (jUUj) and (pTTj) corresponding to our choice 
5 = 0. Thus, we must compose Cj with take a Mellin transform, and extract 

the residue at 0, — i, — 2i in the complex z-plane. To do the Mellin transform (jHOJ), we 
perform first a short distance expansion of the integrand in ()114|) around x 3 , using the 
corresponding expansions ^U] of the quantities a, v n present in H, Ti. This short distance 
expansion will lead to a sum of terms, each of which is a curvature polynomial at x 3 , 
times a Minkowski distribution in £i, £2 and the Riemannian normal coordinates of y. We 
may also set \ — 1; since the residue are independent of the particular choice of Xi as 
proven in Thm. 2. Thus, the computation of the Mellin-transform reduces to ordinary 
Minkowski integrals, times curvature polynomials at x 3 . Only a finite number of these 
terms will give rise to poles at 0, — i, — 2i, and so all others can be discarded. 

Let us consider in detail the pole — 2i. For this pole, the Minkowski space integrals that 
contribute can be reduced to the integral given in jHIj using some standard "iO-identities" 
for distributions: 



C (x 1 ,x 2 ,x 3 ) {1) 



K 1 



2 6 7T 4 p 



Clo ( 2arctan-^-j + Cl 2 f2arctan^-j + Cl 2 f2arctan^- ) + 



, (115) 



3 The usefulness of this operator lies in the identity A(a n ) = 4n(n + l)a r ' 
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where the dots stand for the contributions from the poles 0, — i. Here, we have set 

P= V\^2 + S 2 S 3 + 6 3 6 1 1 (116) 

as well as 

<*i = (6-6)(6-£i), fc = (& -&)(&- 6), *3 = (6-£i)(£i-&), (H7) 

and as before it is assumed that the points & G K 4 are pairwise spacelike, and £ 3 = 0. 
The function Cli is the Clausen function (HOj- The arguments of these functions in 
the above expression are given by twice the angles of a triangle with sides of length 
V (£i — ii) 2 , a/ (£2 — £s) 2 and a/ (£3 — £i) 2 , and p represents the area of that triangle, see 
the figure. 




Thus, the result for Cj manifestly has the simple form claimed in the scaling expansion 
in Theorem 2. There are no curvature terms in the terms that we have displayed, but those 
arise from the other poles 0, — i. These contributions may be obtained in closed form using 
the Minkowski integrals of |Hj. The end result is a sum of terms R^ uap W puap (^i, £2, £3), 
where W£ vap are Lorentz invariant distributions in the Riemannian normal coordinates 
of Xi around x%. However, the expressions for these distributions are rather complicated 
and will therefore be given elsewhere [26 . 

It should be clear that our method is not confined to the above example, but in prin- 
ciple only limited by the ability to perform complicated Minkowski integrals of Feynman 
type. 
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